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We develop a stochastic calculus for processes which are built by convoluting a pure jump, 
zero expectation Levy process with a Volterra-type kernel. This class of processes contains, for 
example, fractional Levy processes as studied by Marquardt [Bernoulli 12 (2006) 1090-1126.] 
The integral which we introduce is a Skorokhod integral. Nonetheless, we avoid the technicalities 
from Malliavin calculus and white noise analysis and give an elementary definition based on 
expectations under change of measure. As a main result, we derive an Ito formula which separates 
the different contributions from the memory due to the convolution and from the jumps. 

Keywords: convoluted Levy process; fractional Levy process; Ito formula; Skorokhod 
integration 

1. Introduction 

In recent years, fractional Brownian motion and other Gaussian processes obtained by 
convolution of an integral kernel with a Brownian motion have been widely studied as a 
noise source with memory effects (see, e.g., Alos et al. [2], Bender [4], Biagini et al. [5] and 
the survey article by Nualart [15]). Potential applications for noise sources with memory 
arc in such diverse fields as telecommunication, hydrology and finance, to mention a few. 

In Marquardt [14], fractional Levy processes were introduced. While capturing memory 
effects in a similar fashion as a fractional Brownian motion does, the convolution with 
a Levy process provides more flexibility concerning the distribution of the noise (e.g., 
heavy tails). In this paper, we consider a larger class of processes by convolution of 
a rather general Volterra-type kernel with a centered pure jump Levy process. These 
convoluted Levy process may have jumps and/or memory effects depending on the choice 
of the kernel. Following the elementary S'-transform approach developed by Bender [4] 
for fractional Brownian motion, we motivate and construct a stochastic integral with 
respect to convoluted Levy processes. The integral is of Skorokhod type and so its zero 
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expectation property makes it a possible choice to model an additive noise. As a main 
result, we derive an Ito formula for these integrals. The Ito formula clarifies the different 
influences of jumps and memory effects, which are captured in different terms. 

The only other paper of which we are aware that treats integration for a similar class 
of processes is [7]. The class of filtered Poisson processes considered in their paper is 
analogously defined by replacing the Levy process by a marked point process in the 
convolution. However, we emphasize that our approach allows the Levy process to be of 
infinite variation and that our Ito formula for the Skorokhod integral is quite different 
from the one of Decreusefond and Savy [7] derive for a Stieltjes integral only. 

The paper is organized as follows. After some preliminaries on Levy processes and con- 
voluted Levy processes in Section 2, we discuss the S'-transform in Section 3. The results 
from Section 3 motivate a definition for a Skorokhod integral with respect to convoluted 
Levy processes which is given in Section 4. In this section, some basic properties of this 
integral are also discussed. Section 5 is devoted to the derivation of the Ito formula, while 
some results are specialized to fractional Levy processes in Section 6. 

2. Preliminaries 

2.1. Basic facts on Levy processes 

We state some elementary properties of Levy processes that will be needed below. For a 
more general treatment and proofs, we refer to Cont and Tankov [6] and Sato [18]. For 
notational convenience, we abbreviate Mq = K \ {0}. Furthermore, ||/|| is the ordinary 
L^-norm of the function / : K — > K and the corresponding inner product is denoted by 
(/iS')l2(r). In this paper, we assume as given an underlying complete probability space 
(17,^, P). Since the distribution of a Levy processes L on {^l,T,P) is infinitely divisible, 
L is determined by its characteristic function in the Levy-Khinchine form E[e™^^*-'] = 
exp{i-0(u)}, t > 0, where 



'4>{u) = l^U - \v^(T^ + / (c' 



1 -iua;l{|^|<i})i^(da;). 



ue R, 



(1) 



R 



7 e M, (T^ > and v \s & Levy measure on M that satisfies 




For any measurable set i? C Mq x let 



N{B)^^{s>Q:{Ls-Ls-,s 



be the Poisson random measure on Rq x R, with intensity measure n(dx,ds) = z^(da;)ds, 
that describes the jumps of L. Furthermore, let N{dx, As) = Nidx, ds) — v{dx) As be the 
compensated jump measure of L (see, e.g., Cont and Tankov [6], Definition 2.18). 
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Assuming that ly satisfies J^^^^i x'^h'{dx) < oo, L has finite mean and variance given by 

var(L(l))= / x^iy{dx)+a^. (2) 

If (7 = in (1), we call L a Levy process without Brownian component. In what follows, 
wc will always assume that the Levy process L has no Brownian part. Furthermore, we 
suppose that E[L(1)] = 0, hence 7 = — ^^^^^-^xvidx). Thus, (1) can be written in the form 

/"(e™^-l-iux)z^(da;), u G M, (3) 

and L = {L(t)}t>o can be represented as 

L{t)^ [ [ xN{dx,ds), teR. (4) 

In this case, L is a martingale. In the sequel, we will work with a two-sided Levy process 
L = {L(t)}tgR, constructed by taking two independent copies {Li{t)}t>o, {i2(i)}t>o of 
a one-sided Levy process and setting 

UiW, iH>0, 
^^^>~\L2i-t-), \ft<0. 

From now on, we will suppose that T is the completion of the cr-algebra generated by 
the two-sided Levy process L and will define L'P{il) := L'p {Q, , T , P) . 

2.2. Convoluted and fractional Levy processes 

We call a stochastic process M = {M(t)}teR given by 

M{t)^ [ /(t,s)L(ds), teR, (6) 
Jr 

a convoluted Levy process with kernel /. Here, / :RxR^R is a measurable function 
satisfying the following properties: 

(i) /(i,-) e L^{R) for aU f e M; 

(ii) f{t, s) = whenever s > t > 0, that is, the kernel is of Volterra type; 
(in) /(O, s) = for almost aU s, hence M(0) = 0. 

Furthermore, we suppose that L = {L{t)}t£R is a Levy process without Brownian com- 
ponent satisfying E[L(1)] = and E[|L(i)|™] < 00 for all rn G N. Hence, the process M 
can be rewritten as 

M{t)^ / / f{t,s)xN{dx,ds), teR. (7) 
Jr Jro 
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Since f(t,-) G ^^(R), the integral (7) exists in L'^{n,P) and 

E[Af(i)^]=E[L(lf] / /2(t,.)d,s = E[L(in||/(i,.)lli2m)- (») 



As an important class of examples of convoluted Levy processes, we now consider 
univariate fractional Levy processes. The name "fractional Levy process" already suggests 
that it can be regarded as a generalization of fractional Brownian motion (FBM). We 
review the definition of a one-dimensional fractional Levy process (FLP). For further 
details on FLP's, see Marquardt [13, 14]. 

Definition 2.1 (Fractional Levy Process (FLP)). Let L = {L{t)}t^]^ be a Levy 
process on M with E[L{1)] = 0, E[L(1)^] < oo and without Brownian component. For 
fractional integration parameter < d < 0.5, a stochastic process 

1 /"^ 

M4t)^^,^^J^Jit-s)i-i-s)i]Lids), teR, (9) 

is called a fractional Levy process (FLP). 
Note that the kernel (9) given by 

satisfies conditions (i)-(iii). Thus, fractional Levy processes are well defined and belong 
to L^{n) for fixed t. 

Moreover, the kernel can be represented by fractional integrals of the indicator function. 
Recall that for < a < 1 , the fractional integral of Ricmann-Liouville type is defined 
by 

r(a) Jx 

(/?/)(x) = -i^ r /(t)(x-tr-Mt 

if the integrals exist for almost all a:; e M. 

In terms of these fractional operators, fractional Levy processes can be rewritten as 

/oo 
{I^X[o,t])(s)L{ds), teR, (11) 
-oo 

where the indicator X[a,6] is given by (a, 6 G R) 

r 1, ifa<t<6, 
X[a,6](t) = < -1, if6<i<a, (12) 
1 0, otherwise. 
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Remark 2.2. The distribution of Md{t) is infinitely divisible for all t <E R, 

Y,[Md{tf] = t2^+iE[i(l)2], t e M, and 

E[exp{izAfdW}] = cxpi / [ {e'^f'^'^-' -l-izft{s)x)v{dx)<la\, t,zGR. (13) 
1,/eJk J 

3. The Levy Wick exponential and the S'-transform 

One of our aims is to introduce a Hitsuda-Skorokhod integral for convoluted Levy pro- 
cesses without touching the technicalities of Malliavin calculus and white noise analysis. 
Our approach is based on the ^-transform, which uniquely determines a square- integrable 
random variable by its expectation under an appropriately rich class of probability mea- 
sures. As a preparation and motivation, we compute the S-transform of Ito integrals with 
respect to the compensated jump measure N in this section. This result then yields a 
simple definition for anticipative integrals with respect to N . 
We begin with some definitions. 

Definition 3.1 (Levy Wick exponential). Let S{M?) denote the Schwartz space of 
rapidly decreasing smooth functions on R^. For r/ G S, where 

7] e S'(R2) : r/(a-, t) > -1, r]{0, t) = 0, ^'?(0, t) = 0, for all t, xe 

the Wiener integral is defined by 

hiv)^ [ [ r^ix,s)N{dx,ds) (14) 

JR JRa 

and the Wick exponential of Ji (rj) by 

exp*(/i(77)) = cxp( / / log{l + iix,t))N{dx,dt)- f f 7?(a:, t)z.(d.T, di) | • (15) 



Remark 3.2. (i) By Theorem 3.1 in Lee and Shih [11] 

/„(r/®") 



expO(/i(,,)) = ^^^i^, (16) 



n 

n=0 



where /„ denotes the multiple Wiener integral of order n with respect to the compensated 
Levy measure. This resprcsentation justifies the name 'Wick exponential'. 

(ii) Since exp*(Ji(?7)) coincides with the Doleans-Dade exponential of Ii{t]) at t = oo, 
it is straightforward that for r;, r) G S, we have 

E[exp<>(/i(77))] = 1 and E[exp*(/i(77)) • exp^(/i(7?))] = exp{(7?,7y)i2(,^^)}, 



where A denotes the Lebesgue measure. 
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We can now define the S'-transform. 

Definition 3.3 (S -transform) . For X G L^{il,P), the S-transform SX of X is an 
integral transform defined on the set S by 



Various definitions of the S'-transform can be found in the hterature, which differ 
according to the chosen subset of deterministic integrands. Our choice of S is particularly 
convenient because of the smoothness of its members. Moreover, it is a sufficiently rich 
set, as demonstrated by the following theorem. It states that every square-integrable 
random variable is uniquely determined by its S'-transform. 

Proposition 3.4. The S-transform is injective, that is, if S{X)(ri) = S{Y)(ri) for all 
rj e S, then X = Y. 

Proof. The assertion is proved in L0kka and Proske [12], Theorem 5.3, by reformulating 
a more general result from Albeverio et al. [1], Theorem 5. □ 

We shall now calculate the S-transform of an Ito integral with respect to the com- 
pensated jump measure N. To this end, let T > and X : Mq x [0, T] x ^ IR be a 
predictable random field (with respect to the filtration !Ft generated by the Levy process 
L{s),0< s<t) satisfying 



The compensated Poisson integral J^^X{y,t)N{dy,dt) then exists in L^{fl,P). 

The following theorem characterizes this integral in terms of the S-transform. The 
result was derived by L0kka and Proske [12], Corollary 7.4, by lengthy calculations in- 
volving multiple Wiener integrals. Here, wc provide a short proof which only makes use 
of classical tools such as the Girsanov theorem. 

Theorem 3.5. Let X denote a predictable random field satisfying the above integrability 
condition. J^^X(y,t)N{dy,dt) is then the unique square-integrable random variable 
with S-transform given by 



(SX)(77)-E'3"[X], 



(17) 



where 



dQ„ = exp*(/i(77))dP. 




■ T 



[ [ SiXiy,t))iriMy,tMdy)dt, rjeE 

Jo JMo 



(18) 
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Proof. Applying Girsanov's theorem for random measures ( Jacod and Shiryaev [8] , The- 
orem 3.17), we obtam that under the measure Q^, the compensator of N(dy,dt) is given 
by {1 + f]{y,t))i'{dy) dt. Hence, 



Xiy,t)Nidy,dt) 



JRo 



X{y,t)ri[y,t)y{dy)dt 



(19) 



is a Q^-local martingale. In particular, if = ti < • • • < < oo is a localizing sequence 
of stopping times with lim^y^cx) Tjy = oo a.s., then 



lim E'5i 

N^oo 



TArjv 



Xit,y)Nidy,dt) 



hm 



TAr„ 



X{t,y)r,{y,t),^{dy)dt 



Xit,y)rjiy,tyidy)dt 



by a straightforward application of the dominated convergence theorem. 
To treat the limit in the first line, note that 



TAr„ 



= E^ 



cxpO(/i(77)) 



TAt„ 



X{y,t)Nidy,dt) 



Xiy,t)Nidy,dt) 
The integrand on the right-hand side is dominated by 

exp^(/i(77)) sup / / X{y,t)N{dy,dt) 

Q<u<T Jo jRn 



which is P-integrable by Holder's inequality, Doob's inequality and the assumed integra- 
bility of the random field. Thus, 



Xiy,t)N{dy,dt) 



JRo 



X{y,t)r^{y,tMdy)dt 



JRo 



and the assertion follows by applying Fubini's theorem. 

Note that the last identity shows that the (5,,-local martingale (19) is indeed a 
Q^-martingale. □ 

Example 3.6. By separately applying the previous theorem to both sides of the two- 
sided Levy process, we derive 



SiMit)M 



f{t,s)yri{y,s)v{dy) ds 
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since 

M{t)^ f f{t,s)L{ds)= f I f{t,s)yN{dyAs). 

The S'-transform characterization in the previous theorem gives rise to a straightfor- 
ward extension to anticipative random fields. 

Definition 3. 7. Suppose X is a random field. 

(i) The Hitsuda-Skorokhod integral of X with respect to the compensated jump mea- 
sure N is said to exist in L^(f2) if there is a random variable <f> € L^(ri) such that for all 

i-T 

, . Jdt. 

/O JRo 



S^{vi)^ r [ S{X{y,t)Mrj{y,tMdy), 

Jo JRn 



It is denoted by <^ = J^^X{y,t)N^{dy,dt). 

(ii) The Hitsuda-Skorokhod integral of X with respect to the jump measure N is 
defined as 

X{y,t)N\dy,dt).= j [ X{y,t)N\dy,dt)+ f [ Xiy,tHdy)dt 
JRo Jo Jro Jo Jro 

if both integrals on the right-hand side exist in L'^{fl). 

Remark 3.8. From the previous definition, we get immediately that 

s{j I X{y,t)N%dy,dt))irj)^ [ [ SiXiy,t))i7j){l + Tjiy,t)Hdy)dt. 
\Jo JRo / "'O JMo 

Clearly, if the integrand is predictable, this Skorokhod integral reduces to the ordinary 
stochastic integral for random measures and the diamond can be omitted in this case. 

Remark 3.9. Theorem 3.5 implies that 

ft 



S{N{A,[OMi.v)= I I V{v.s)y{dy) 

Jo J A 



ds. 



Hence, we can write, in a suggestive notation, 

s(f( X{y,t)N'>{dy,dt)){ij)= r f S{X {y,t)MS{N{dy,dt)){^). 
\Jo JRo ' Jo JRo 

In view of Example 3.6, Theorem 3.5 can be specialized to integrals with respect to 
the Levy process L as follows. 
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Corollary 3.10. Let < a <b and X : [a, b] x Q ^ M. be a predictable process such that 
E[j'' dt] <oo. X{s)L{ds) is then the unique square-integrable random variable 
with S-transform given by 

' f SiXmv)-^S{Lmv)dt, ryes. 

We close this section with a remark concerning the relationship between the Skorokhod 
integral with respect to N and ordinary integration. 

Remark 3.11. Suppose u{y,t) is a simple random field of the form 

=-FlAx(a,f,](2/,i), a<beR,FeL\n), 
where A C Rq satisfies iy{A) < oo. Then, in the sense of ordinary integration, 

uiy,t)Nidy,dt)=FNiA, (a, 6]). 



We want to relate this expression to Skorokhod integration with respect to N. We shall 
suppose that F = exp* (/(/)) for some f E L^{i^x X), where, in generalization of Definition 
3.1, 

exp^(/(/)):=exp( / / f{y,t)N{dy,dt)\l[{l + f{AL{s),s))e'f^^^^-^^'^\ 

A direct calculation then shows that ior f E L^{i' x X), r/ E S, 

expO(/(/)) exp^(/(7?)) = exp^(/(/ + r/ + frj)) cxpj^ ^ f{y, tHy, tHdy) dt 

Consequently, by a slight extension of Remark 3.9, 

s( [ [ u{y,t)N{dy,dt)]{Tj) 
\JrJro ' 

= exp(/ / /(y,t)77(2/,t)Kdy)dAE[expO(/(/ + r; + /,7))iV(A,(a,6])] 



(5exp-(/(/)))(77) /'/ [7?(y,s) + /(y,s) + r;(y,s)/(y,s)]:.(dy)ds 

Ja J A 



Si uiy,t)fiy,tMdy)dt (tj) 

\JrJro ' 

+ s( [ [ [u{y,t)f{y,t) + u{y,t)]N'>{dy,dt)\{rj). 
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Let us now define the Malliavin derivative of a Wick exponential by 

Z?,.,exp^(/(/)) = /(y,s)exp^(/(/)), 

wliich can be extended to a linear closed operator acting on a larger class of random 
variables (see, e.g., Nualart and Vivos [16]). We then arrive at the formula 



u{y,t)N{AyAt)^ / / Dy^tu{y,t)v{dy)As 

[u[y,t) + Dy^tu{y,t)]N\Ay,ds). 



We conjecture that this formula can be extended by approximation to a larger class of 
random fields. 



4. A Skorokhod integral for convoluted Levy 
processes 

In this section, we define the Skorokhod integral for convoluted Levy processes and state 
some basic properties. The definition is strongly motivated by Corollary 3.10 above. 

Definition Suppose that the mapping 

t^SiMit)M 

is differentiate for every rj Suppose B (iM. is a Borel set and X : B x ^M. is a 
measurable stochastic process such that X{t) is square-integrable for each t Cz B . X is 
said then to have a Hitsuda-Skorokhod integral with respect to M if 

5(X(.))(r/)^5(Af e L\B) for anyr,GE 
and there is a $ G L'^(n) such that for all 7? G 

smv) = l^s{x{t)M^^s{M{t)Mdt. 

In that case, $ is uniquely determined by the injcctivity of the 5'-transform and we 
write 

$ = / X{t)M''{dt). 

Remark 4-2. (i) Lemma 5.1(ii) below provides some sufficient conditions for the differ- 
entiability of the mapping f i— > S{M{t)){if} in terms of the convolution kernel. 
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(ii) The definition of the Skorokhod integral does not require conditions such as pre- 
dictabihty or progressive measurabihty. Hence, it also generalizes the ltd integral with 
respect to the underlying Levy process to anticipative integrands. 

(iii) Since the Levy process itself is stochastically continuous, the S'-transform cannot 
distinguish between L{t) and L{t—) for fixed t. Consequently, we obtain, for example, 

f L(s)i*(ds)= / L(s-)i*(ds)= / L(s-)i(ds), 
Jo Jo Jo 

where the last integral is the classical Ito integral. 

The following properties of the Skorokhod integral are an obvious consequence of the 
definition. 

Proposition 4.3. (i) For alla<beE., M (b) - M (a) = jj^ APidt). 

(ii) Let X : B X fl ^ L^(r2) be Skorokhod integrable. Then 

X(t)Ar(dt)= / lBit)X{t)M''{dt), 
Jr 

where 1b denotes the indicator function of the set B. 

(iii) Let X:Bxn-y L'^ifl) be Skorokhod integrable. Then E[/g X(t)Af^(dt)] = 0. 

We note that (iii) holds since the expectation coincides with the S'-transform at = 0. 
The zero expectation property makes the integral a promising candidate for modeling an 
additive noise. 

Example 4-4- As an example, we show how to calculate AI{t)M^{dt). In the fol- 
lowing manipulations, N'^ denotes the compensated jump measure under the probability 
measure Qn = cxp*(/i (77)) dP. In particular, it follows from Girsanov's theorem, as in 
the proof of Theorem 3.5, that 

M{T)^ f [ f{T,s)yN''{dy,ds)+ [ [ f{T,s)yT^{jj,s)u{dy)ds. 

J-ocJs.0 J-ooJRo 

By this identity, integration by parts and Example 3.6, we obtain 
5(^2^'^M(i)Ar(dt)^(r?) 

^2 1^ SiMmv)^^S{M{t)Mdt 

iS{M{T))ir,)f=[ I I fiT,s)yrjiy,sMdy)ds' 
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'' f{T,s)yN\dy,ds)+ f f f{T, s)yr,iy, sHdy) ds) 



— OO J Kq J —oo 

/ f{T,s)yN^{dy,ds 

-oo JRo 



= S{M{TY){ij)- I / /(T,s)V(l + ^(y,s))Kd2/)ds 



oo JKa 



Here, we have used the fact that /Tqo/r f{T,s)yN^{dy,ds) has zero expectation and 



variance f{T, s)^y^(l + 77(1/, s))i/(dy) ds since the compensator of N under is 

given by (1 + ?/(?;, s))v{dy, ds). 

Hence, from Remark 3.8, we derive the identity 



T pT 

2 



2/ M{t)M^{dt)=M{TY~ / /(T,s)VA^(dy,ds) 

Jo J-ooJKo 



-oo<s<T 

provided both members on the right-hand side exist in L'^iVl). 

Remark 4-5. Applying the same techniques as in Example 4.4, one can easily obtain, 
for a <b, 

M{a) flM%dt)^ f M{a)M%dt)+ T / f{a,s){f{h,s) - f{a,s))y^N{dy,ds). 

Ja Ja Jo JMo 

Hence, ordinary multiplication with a random variable, which is measurable with respect 
to the information up to the lower integration bound, cannot, in general, be introduced 
under the integral sign if the kernel depends on the past. 

5. Ito's formula 

In this section, we will derive an Ito formula for convoluted Levy processes. The proof 
is based on a calculation of the time derivative of S{G{M{t))){ri). It may be seen as a 
generalization of the calculations in Example 4.4. This technique of proof is in the spirit of 
Kubo [9], Bender [3] and Lec and Shih [10], where this approach was applied to obtain Ito 
formulas for generalized functionals of a Brownian motion a fractional Brownian motion, 
and a Levy process with Brownian component, respectively. 

During the derivation of the Ito formula, we have to interchange differentiation and in- 
tegration several times. Under the following (rather strong) conditions on the convolution 
kernel, these manipulations are easily justified. However, the Ito formulas below may also 
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be viewed as generic results which hold for more general kernels (with the technicalities 
to be checked on a case-by-case basis). 
We make the following assumptions: 

(HI) there are constants a < < 6 such that supp(/) C [a, 5]^; 

(H2) / is continuous and bounded on [a, 6]^ \ {{t, s); i = s}; 

(H3) linisjf /(t, s) = f(t, t) and the mapping 1 1-^ /(t, t) is continuous; 

(H4) / is continuously differentiablc on (a, b)^ \ {{t, s);t ^ s} with bounded derivative. 

Lemma 5.1. Under (H1)-(H4), we have the following: 

(i) Fora<t<b, 

M{t) = f{t,t)L{t) - f{t,a)L{a) - ^ L{s)^J{t,s)ds. (20) 

In particular, M(t) has a modification which is RCLL and stochastically continuous. 
Moreover, 

AM{t) = f{t,t)ALit). 

Hence, M is continuous on [a, b] if and only if f{t, t) ~ for all a <t <b. 

(ii) The mapping [a, 6] ^ M, 1 1-^ {SM{t)){r]) is continuously differentiablc for all rj 
and 

-^(5M(t))(77)= r / ^f{t,s)yf^{y,sHdy)ds + f{t,t) f yr^{y,tMdy). (21) 

J-OD J Mo CIS JRo 

Proof, (i) Fix a modification of L which is right-continuous with left limits (RCLL). 
Formula (20) follows from the definition of M and integration by parts, which is justified 
by (H3)-(H4). The second and third terms on the right-hand side are continuous in t 
by (H2) and the boundedness of -^fit^s), respectively. The first term is stochastically 
continuous and RCLL since L has these properties and ti—ff(t,t) is continuous. The 
other assertions in (i) are obvious consequences. 

(ii) can easily be obtained by differentiating the expression in Example 3.6. □ 

Example 5.2. The following prominent examples satisfy conditions (H1)-(H4): 

1. one-sided shot noise processes defined by the kernel 

^ ^ ' ' \ 0, otherwise, 

for constants T* > and fc; 

2. one-sided Ornstein-Uhlenbeck type processes defined by the kernel 



^~k(t-s}^ 0<s<t<T* 
0, otherwise, 
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for constants T* > and fc > 0. 

From the previous lemma, one directly obtains that the shot noise processes have con- 
tinuous paths, while the Ornstein-Uhlenbeck-type processes exhibit jumps. 



as 



To state the Ito formula precisely, we finally recall that the Wiener algebra is defined 

{GeL\R);TGeL\R)}, 

where T denotes the Fourier transform. Note that the space of rapidly decreasing smooth 
functions is included in the Wiener algebra. 

The first version of Ito's formula requires that the underlying Levy process is a finite 
variation process. 

Theorem 5.3 (Ito formula I). Let (H1)-(H4) hold, 0<T<b and 

|x|j/(da;) < oo. 



Furthermore, assume that G £ C"'^(]R) with G, G' G A{M.) bounded. Then 



T / i-t 







G'(M{t) + xf{t, s))x^fit, s)N''{dx, ds) ] dt 



exists in L'^{fl) and 

G{M{T)) = G(0) - (^^ x,y{dx)^ £ G"(M(i)) (^/(t, t) + j' ^J{t, s) ds) dt 
+ ^ G{M{t))-G{M{t-)) 



0<t<T 
i-T / rt 







G'{M{t-) + xf{t, s))x^f{t, s)N''{dx, ds) ] dt. 



In the general case, the Ito formula reads as follows. Indeed, the previous formula can 
be derived from the general one by rearranging some terms. 

Theorem 5.4 (Ito formula II). Let (H1)-(H4) hold and 0<T<b. Furthermore, 
assume that G G C^(M) with G, G' € A(R). Then 

G{M{T)) 

= G(0)+ / G'{M{t~))Ae{dt) 
Jo 
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+ G{M{t))^G{M{t-))~G'{M{t-))AM{t) 

0<t<T 

+ (^f £ {G'{M{t-) + s)) - G'{M{t-)))x^J{t, s)iV^(dx, ds) j At, 

provided all terms exist in Lp'iVl). 

We would like to emphasize that the Skorokhod integrals with respect to N in the above 
versions of Ito's formula do not, in general, reduce to ordinary integrals for the following 
reason. The time variable of the Skorokhod integral is s, but the integrand depends 
on M through the value M{t~), where t> s. Therefore, the integrand is typically not 
predictable as a process in the variable s. 

The above versions of Ito's formula (but not their exact assumptions) reduce to well- 
known formulas for Levy processes with the choice /(t, s) = X(o,t](s) as, in this case, the 
last Skorokhod integral with respect to N vanishes. We recall that M has independent 
increments if and only if ^f{t,s) ~ for all t. Hence, the contributions from disconti- 
nuities and memory effects are nicely separated in the above Ito formulas. Finally, note 
that the formula for M(tY from Example 4.4 can be recovered by formally applying the 
Ito formula II with G{y) ~ y^. 

Remark 5.5. ltd formula II has the drawback that the conditions do not guarantee 
that all members of the identity exist in L'^{n). However, the manipulations below can 
be recast in a white noise framework, as developed in [17], in a way that all members 
exist as generalized random variables. 

The remainder of this section is devoted to the proof of the Ito formulas. As a gen- 
eral strategy, we wish to show that both sides of the asserted identities have the same 
^-transform. Indeed, the following calculations show how to identify the right-hand side 
constructively. We first write 

S{G{M{T))M = G(0) + ±SiGiMit))M dt 

and then calculate ■j^S{G{M{t))) explicitly. To achieve this, we apply the inverse Fourier 
theorem and obtain, for G S ^(K), 

S{G{M{t))){T]) = E'^" [G[M{t)]] = I .FG(u)E'3" [e™^^(*)] du. (22) 

To differentiate this expression, we calculate the characteristic function of M under Q^. 

Proposition 5.6. Let AI ~ {M{t)}t£R be a convoluted Levy process as defined in (6), 
with kernel function f. Then 

5(e'"*^(*))(7?) 
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explains {M{t)){T]) 



1 — mxf{t, s))(l + s))]iy{dx) ds 



Proof. It follows from the proof of Theorem 3.5 that 



L^{t):^L{t)- I / x-q{x,s)v{dx)ds 



is a (5,;-martingale with zero mean. Applying Girsanov's theorem for semimartingales 
(Jacod and Shiryayev [8], Theorem 3.7) yields that L'^ has semimartingale characteris- 
tics (7^,0, v'^), where 7? = — x{\ + i^{x, s))z^(da;) and ^^{dx) = (1 + rj{x, s))v{dx). 
Hence, 



5(exp{itiL'3(0})(?7) = 



exp 



1 — iux][l + rj[x, s)]v[dx) ds > . 



JRo 



Finally, 



5(exp{iuAf(i)})(77) 



exp'^ iit / f{t,s)L(ds) 
J —00 

iu / f{t,s)L^{ds)+m / f{t,s) / xr]{x,s)v{dx)ds 



exp 



: exp 



X exp 



00 JKq 

t 



^f(t,s) -i-.[uxf{t,s)][\+r]{x,s)]v{dx) ds 



luxfit, s)rj{x^ s)v{dx) ds 



Taking the S'-transform of M into account, which was calculated in Example 3.6, the 
assertion follows. □ 

By introducing the derivative under the integral sign, we get 
d 



dt 



^gQ,[gi«A/(t)] f [(e™-/(*-t)_i_iy2./(i^t))(l + ?7(a;,f))]i/(dx) 



i„^_/(t,5)(e'"-^(*'^)-l)(l + 77(x,<)) 



v{dx) ds 



4_E'3''[e'"*^W]iw^^(M(i))(7?). 



(23) 



Stochastic calculus for convoluted Levy processes 515 

Combining (22) with (23) and again interchanging differentiation and integration 
(which can be justified under (H1)-(H4) since G,G' £ A{M.)), we obtain 

^SiGiAimiv) 

^ / (^G)(M)EQ"[e'"^^(*)] / [(e™"-^(*'*)-l-iux/(i, <))(! + 7y(a;,t))];.(da;)du 



^ / (J^G)(w)E'3'>[e'"*^(*)] 
2n Jr 



iy^A/(i,s)(e'"-^(*'^)-l)(l+r?(x,t)) 

+ ^ / (^G)(w)E'?''[e'"^^W]iu^5(M(i))(7?)dM 
V27T 



i/(da;) ds du 



= :(/) + (//) + (///). 
Standard manipulations of the Fourier transform, together with (22), now yield 

(/) = -^/ / [{J^G{- + xf{t,t)){u)^{J^G){iL)-xf{t,t){J^G'){u)] 



27T 

X E'^'^ [e™^^(*)] (1 + 7y(a;, t))iy{dx) du 
S{G{M{t~)+xf{t,t)) - G{M{t~)) - xf{t,t)G'{M{t^))){ij) 

X (1 + ri{x, t))iy{dx). 
The second term can be treated analogously and thus, 

i")^ f I x^f{t,s)S{G'{M{t~) + xf{t,s))^G'iM(t-mrj) 
X (1 + ri{x, t))iy{dx) ds. 

Finally, (///) = 5(G'(M(t-)))(77)^5(M(t))(77). 

We now collect terms and integrate t from to T, whence 

5(G(Af(T)))(77)-G(0) 

/ S{GiM{t~) + xf{t, t)) - G{M{t-)) - x/(i, t)G\M{t-m'n) 

JRo 

X {l + r^{x,t))u{dx)dt 



x-f{t, s)S{G'{Mit-) + xf{t, s)) - G'(M(t-)))(7y) 
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X {l + r]{x,t))iy{dx)dsdt 

+ 1^ 5(G'(Af(i-)))(r;)^5(M(i)) W dt 
=:{i) + {U) + iui). (24) 

From Remark 3.8, we get 

{i) = s( r [ G{M{t~) + xf{t, t)) - G{M{t-)) - xf{t, t)G'{M{t~))N'>{dx, dt)) (t?) 
\Ja JRo / 

= s( J2 G{M{t))-G{M{t-))-G'{M{t-))AAI{t)]{rj), 



KO<t<T 



where the second identity holds because the Skorokhod integral is an Ito integral by 
predictability (and by taking Lemma 5.1(i) into account). Similarly, 

f ^ x^J{t,s)[G'{M{t-)+xf{t,s))-G\M{t~))]N%dxM<^t){l)- 
Finally, by the definition of the Skorokhod integral with respect to M, 



{tii) = sl^j^ G'{M{t-))M\dt)yri). 



Hence, both sides of Ito formula II have the same S-transform, which proves this formula. 
To get Ito formula I, we rearrange the terms in (24). By Lemma 5.1(ii), 



d „ „, . /• /•* d 



S{M{t)){r]) ^ f{t,t) xr]{x,t)v{dx) + 3-/(^,5) / xri{x, s)iy{dx) ds 



./Ro ./-nn dt 



Thus, by (24) and similar considerations as above. 



7* 



x^ fit, s)S{G'{M{t-) + xf{t, s)))(v)(l + v{x, t)),y{dx) dsdt 



s 



G{M{T)) - G(0) + (^^ x,y{dx)) £ G'{M{t)) (^f{t, t) + j' ^^f{t, s) ds^ dt 
- G(M(t-) + AM(t))-G(Af(t-)))(,7). 



o<t<r 
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The expression under the ^'-transform on the right-hand side clearly belongs to L^(ri) 
under the assumptions of Ito formula I. Then, by Remark 3.8, the Skorokhod integral 

^ f [ x^f{t,s)G'{M{t-)+xf{t,s))N'>{dx,ds)dt 

J-ocJs.0 

exists in L'^{fl) and coincides with the expression under the 5-transform on the right- 
hand side. This proves Ito formula I. 



6. Stochastic calculus for fractional Levy processes 

We shall now specialize from a convoluted Levy process to a fractional one. In Marquardt 
[14], a Wiener- type integral with respect to a fractional Levy process is defined for de- 
terministic integrands. Its domain is the space of functions g such that 5 e L^(M) and 
it can be characterized by the property 



Jm. 



g{s)MMs)= / {llg){s)Lids). 



The following theorem shows that a similar characterization holds for Skorokhod integrals 
with respect to fractional Levy processes. Hence, it also proves, as a by-product that the 
Wiener-type integral is a special case of the Skorokhod integral. 
In preparation, note that 

S{Mdmv)= f f I^X[o,t]{^)yv{s,yHdy)ds. 

JR JRo 

Hence, by Fubini's theorem and fractional integration by parts, we obtain the following 
theorem. 

Theorem 6.1. Suppose Md is a fractional Levy process with < d < 0.5. Then, for all 
±SiAUt)){ri)= I {Iiil){t,y)yv{dy), 

J Mo 

where, by convention, fractional integral operators are applied only to the time variable 
t. 

Furthermore, suppose that X e LP {R, L'^ (Q)) with p ^ {d + 1/2)-^ . Then 

X{t)Amt)= I {liX){t)L%dt) 

in the usual sense, that is, if one of the integrals exists, then so does the other and both 
coincide. 
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Proof. The proof follows the same lines as that of Theorem 3.4 in Bender [4]. □ 

Note that only Ito formula II makes sense for fractional Levy processes. When we 
formally apply this Ito formula, the following observation is noteworthy. For d> 0, the 
process Af^; is continuous and has memory, whence 

G{Md{T)) ^ GiO) + f G'(Mrf(t-))Ar(dO 



- G'(Mrf(t-)) j ^(t - sf-'N^idx, ds)j dt. 

However, the Levy process L itself comes up as limit of Md when d tends to 0. As this 
process has independent increments and jumps, its well-known Ito formula reads 

G(L(r)) = G(0)+ r G'{L{t-))L{dt) 



J2 Gim) G{L{t-)) G'{L{t-))M{t). 



0<t<T 

So, apparently, the Ito formulas do not transform continuously into each other when 
passing to this limit. This is in sharp contrast to the Gaussian case, in which the Ito 
formula for Brownian motion is recovered by substituting H — 1/2 (the Hurst parameter 
corresponding to d via d = H — 1/2) into the Ito formula for fractional Brownian motions 
(see, e.g., [4]). 



Acknowledgements 

The paper benefited from the constructive comments of two anonymous referees. In 
particular, their remarks helped to clarify the conditions required for the proof of Ito's 
formula. 



References 

[1] Albeverio, S., Daletsky, Y., Kondratiev, Y.G. and Streit, L. (1996). Non-Gaussian infinite 
dimensional analysis. J. Fund. Anal. 138 311-350. MR1395961 

[2] Alos, E., Mazet, O. and Nualart, D. (2001). Stochastic calculus with respect to Gaussian 
processes. Ann. Probab. 29 766-801. MR1849177 

[3] Bender, C. (2003). An Ito formula for generalized functionals of a fractional Brownian mo- 
tion with arbitrary Hurst parameter. Stochastic Process. Appl. 104 81-106. MR1956473 



Stochastic calculus for convoluted Levy processes 



519 



[4] Bender, C. (2003). An S-transform approach to integration with respect to a fractional 

Brownian motion. Bernoulli 9 955-983. MR2046814 
[5] Biagini, F., 0ksendal, B., Sulem, A. and Wallner, N. (2004). An introduction to the white 

noise theory and Malhavin calculus for fractional Brownian motion. Proc. R. Soc. Land. 

Ser. A Math. Phys. Eng. Sci. 460 347-372. MR2052267 
[6] Cont, R. and Tankov, P. (2004). Financial Modelling with Jump Processes. Boca Raton, 

FL: Chapman and Hall. MR2042661 
[7] Decreusefond, L. and Savy, N. (2006). Anticipative calculus for filtered Poisson processes. 

Ann. Inst. H. Pomcare Probab. Statist. 42 343-372. MR2219714 
[8] Jacod, J. and Shiryaev, A.N. (2003). Limit Theorems for Stochastic Processes, 2nd ed. 

Berlin: Springer. MR1943877 
[9] Kubo, I. (1983). Ito formula for generalized Brownian functionals. Lecture Notes in Control 

and Inform. Sci. 49 (G. Kallianpur, ed.) 155-166. Berlin: Springer. MR0799940 
[10] Lee, Y. and Shih, H. (2000). Ito formula for generalized Levy functionals. In Quantum 

Information II (T. Hida and K. Saito, eds.) 87-105. River Edge, NJ: World Scientific. 

MR1880131 

[11] Lee, Y. and Shih, H. (2004). The product formula for multiple Levy-Ito integrals. Bull. 
Inst. Math. Acad. Sin. (N.S.) 32 71-95. MR2069183 

[12] L0kka, A. and Proske, F. (2006). Infinite dimensional analysis for pure jump Levy processes 
on the Poisson space. Math. Scand. 98 237-261. MR2243705 

[13] Marquardt, T. (2006). Fractional Levy processes, CARMA processes and related topics. 
Ph.D. thesis, Munich Univ. Technology. 

[14] Marquardt, T. (2006). Fractional Levy processes with an application to long memory mov- 
ing average processes. Bernoulli 12 1090-1126. MR2274856 

[15] Nualart, D. (2003). Stochastic calculus with respect to the fractional Brownian motion and 
applications. Contemp. Math. 336 3-39. MR2037156 

[16] Nualart, D. and Vives, J. (1990). Anticipative calculus for the Poisson process based on the 
Fock space. Seminaire de Probabilites XXIV 154-165. Lecture Notes in Math. 1426. 
Berhn: Springer. MR1071538 

[17] 0ksendal, B. and Proske, F. (2004). White noise for Poisson random measures. Potential 
Anal. 21 375-403. MR2081145 

[18] Sato, K. (1999). Levy Processes and Infinitely Divisible Distributions. Cambridge: Cam- 
bridge Univ. Press. MR1739520 

Received March 2007 and revised October 2007 



